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We study the surface scaling behavior of a semi-infinite <i-dimensional O(N) spin system in the presence of 
quenched random field and random anisotropy disorders. It is known that above the lower critical dimension 
die = 4 the infinite models undergo a paramagnetic-ferromagnetic transition for N > N c (N c = 2.835 for 
random field and N c — 9.441 for random anisotropy). For N < N c and d < di c there exists a quasi-long- 
range ordered phase with zero order parameter and a power-law decay of spin correlations. Using functional 
renormalization group we derive the surface scaling laws which describe the ordinary surface transition for 
d > di c and the long-range behavior of spin correlations near the surface in the quasi-long-range ordered phase 
for d < d\ c . The corresponding surface exponents are calculated to one-loop order. The obtained results can be 
applied to the surface scaling of periodic elastic systems in disordered media and amorphous magnets. 

PACS numbers: 68.35.Rh, 64.60.ae, 75.30.Kz 
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I. INTRODUCTION 



The phase diagram and critical properties of spin systems 
with quenched disorder attracted considerable interest for 
decades. One usually distinguishes two types of quenched 
disorder: (i) random-temperature like disorder correspond- 
ing to randomness coupled to the local energy density as, for 
example, in diluted ferromagnets [1]; (ii) random field like 
disorder corresponding to the case when the order parame- 
ter couples to a random symmetry breaking field [2]. The 
influence of random-temperature disorder is rather well un- 
derstood. There exist several powerful methods to study the 
phase behavior and criticality such as perturbative renormal- 
ization group (RG). The effect of the random field disorder 
being more profound is much less studied. The prominent ex- 
ample is the critical behavior of the random field Ising model 
(RFIM) which complete understanding is still lacking despite 
significant numerical, analytical and experimental efforts [3]. 
It has been found that the perturbative calculations including 
standard RG methods lead to incorrect results, in particular, 
to the so-called dimensional reduction (DR). Analysis of the 
Feynman diagrams giving the leading singularities (|4j] or us- 
ing supersymmetry Q predicts that the critical behavior of the 
RFIM in d dimension is the same as that of the pure system in 
d — 2 dimensions. Consequently, the lower critical dimension 
of the RFIM below which there is no true long-range order is 
expected to be = 3. However, the simple Imry-Ma argu- 
ments show that the lower critical dimension of the RFIM is 
in fact di c = 2 [2]. The deviation from the DR prediction is 
also confirmed by the high-temperature expansion [6] and real 
space RG |7|]. The failure of DR can be explained by compli- 
cated energy landscape which renders the perturbation theory 
spoiled to all orders by unphysical averaging over multiple 
minima and maxima. The latter can be formulated in terms 
of supersymmetry or replica symmetry breaking JH 0. To 
overcome this obstacle one needs a non-perturbative method 
or correct resumming the perturbation theory. 

The considerable progress has been achieved last years in 
studying the O(N) models in which disorder couples to the 
iV-component order parameter either linearly as in the ran- 
dom field (RF) case or bilinearly as in a random anisotropy 



(RA) system. These models are relevant for diverse physical 
systems including amorphous magnets 1 10], diluted antifer- 
romagnets in a uniform external magnetic field ifTHl . liquid 
crystals in porous media 012l Il3p . nematic elastomers Q14ll . 
critical fluids in aerogels IU5I417II . vortices in type II super- 
conductors OH], and stochastic inflation in cosmolo gy fl9ll . 
Similar to the RFIM these models suffer of DR |H 12011. It 
was shown that the expansion around the lower critical di- 
mension of the the RF O(N) model d\ c — 4 generates an infi- 
nite number of relevant operators which can be parameterized 
by a single function [21]. However, the RG flow of this func- 
tion has no analytic fixed point (FP). Only almost two decades 
later, being inspired by the progress in disordered elastic sys- 
tems Il22l - l25ll . it was realized that the scaling properties of the 
RF and RA systems are encoded in a nonanalytic FP l26ll . 
The non-analytic FPs control the paramagnetic-ferromagnetic 
phase transitions in the RF and RA O(N) model and allow 
one to compute the critical exponents within e = d — 4 ex- 
pansion II27I1 . The obtained exponents are different from the 
DR prediction. The FRG calculations have been extended to 
two-loop order ||28[ l29ll and the effect of long-range disorder 
correlations has been studied lfl7l [30h . Using developed in 
Ref. 13111 truncated exact FRG it was argued that spontaneous 
breaking of the supersymmetry which leads to a breakdown of 
DR, occurs only below a critical dimension doR ~ 5.1 |32]. 

A more peculiar issue concerns the phase diagram of the 
RF and RA models below d\ c . It is known that for the 
RF model and models with isotropic distributions of random 
anisotropics true long-range order is forbidden below d\ c — 
4 (for anisotropic distributions, long-range order can occur 
even below di c |33]). Nevertheless, quasi-long-range order 
(QLRO) with zero order parameter and an infinite correlation 
length can persist even for d* c (N) < d < d\ c , where d* c (N) is 
the lower critical dimension for the paramagnetic-QLRO tran- 
sition. For example, the Gaussian variational approximation 
predicts that the vortex lattice in disordered type-II supercon- 
ductors can form the so-called Bragg glass exhibiting slow 
logarithmic growth of displacements Q4I1 . This system can be 
mapped onto the three dimensional RF 0(2) model, in which 
the Bragg glass corresponds to the QLRO phase. Indeed, for 
N < N c and d < d\ c , the FRG equations have attractive FPs 
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which describe the QLRO phases of RF and RA models 112611 . 
Despite that the question of the lower critical dimension of 
the paramagnetic-QLRO transition is still controversial. In 
order to study the transition between the QLRO phase and the 
disordered phase using FRG, one has to go beyond the one- 
loop approximation. The truncated exact FRG rf3lll and the 
two-loop FRG [29] performed using a double expansion in 
■v/jej and N — N c provide an additional singly unstable FP 
which controls the transition. Both methods give qualitatively 
similar pictures of the FRG flows: the critical and attractive 
FPs merge in some dimension d^ c (N) < d\ c which is con- 
sidered as the lower critical dimension of the paramagnetic- 
QLRO transition. For the RF 0(2) model, both methods give 
approximately the same estimation d* c ps 3.8(1), and thus, 
suggest that there is no Bragg glass phase in d = 3. However, 
one has to take caution when extrapolating results obtained for 
small \ f\e\ and N — N c . Moreover, in contrast to the model 
of Refs. B26I1 and 12911 which belongs to the so-called "hard- 
spin" models, the system studied in Ref. 113 111 corresponds to 
"soft spins". They can belong to different universality classes 
since the soft spin model allows for topological defects which 
destroy the QLRO. 

The real systems, usually, are finite and have boundaries 
which effect is twofold: (i) the free energy of the system in 
addition to the bulk contribution proportional to the volume 
acquires a new term proportional to the area of the surface; 
(ii) the presence of boundaries breaks the translational invari- 
ance. In general, this can modify the behavior in the boundary 
region extended in the bulk only over distances of the order 
of the bulk correlation length. However, at the bulk critical 
point or in the QLRO phase, the bulk correlation length is in- 
finite so that one can expect that the effect of boundaries to be 
more pronounced. Indeed, the presence of the boundaries in- 
troduces a whole set of critical exponents describing the scal- 
ing behavior at and close to the boundary at criticality ||35ll . 
Several different classes of the surface transitions are known 
depending upon boundary conditions ll36ll . The ordinary tran- 
sition corresponds to the case when the surface magnetization 
is suppressed due to reduced number of close neighbors near 
the boundary so that the surface ordering is completely driven 
by the bulk magnetization. If for some reason the coupling 
between spins on the surface is sufficiently enhanced with re- 
spect to the bulk coupling or there is an external surface mag- 
netic field, the surface may order before the bulk does. The 
latter is called the surface transition. Then the system can un- 
dergo the so-called extraordinary transition in the presence of 
ordered surface. The two lines of the extraordinary transition 
and the surface transition meet at the multicritical point which 
is called the special transition. The last three transition can 
take place only if the dimension of the surface d — 1 is above 
the lower critical dimension for the transition. These transi- 
tions have been studied for various systems with discrete and 
continuous symmetries using different methods, such as RG 
and numerical simulations (for review see ll3~5ll3~7[|38ll '). 

The effect of weak random temperature like disorder on 
the surface criticality was studied using RG methods in 
Refs. i39l l40ll . However, not so much is known about the 
surface criticality in systems with RF disorder. The phase di- 



agram of the 3D semi-infinite RFIM as a function of the ratio 
of bulk and surface interactions and the ratio of bulk and sur- 
face fields has been studied using a mean field approximation 
in Ref. 14111 . The surface criticality of the RFIM has been 
studied numerically in Ref. H42I1 . It was also shown that the 
RF disorder on the surface of a 3D spin system with continu- 
ous symmetry destroys the long-range order in the bulk, and, 
instead, a QLRO emerges |43J. In this work we address the 
question how do the RF and RA disorder in the bulk effect the 
behavior of spin systems with continuous symmetry in vicin- 
ity of free surfaces. In particular, we consider the ordinary 
surface transition of the RF and RA O(N) models for d > 4 
and the spin correlations in the QLRO phase near a free sur- 
face for d < 4. 

The paper is organized as follows. Section|II]introduces the 
model. In Sec.|III|we renormalize the theory and derive the 
scaling laws. In Sec. ||V]we calculate the surface critical ex- 
ponents to one-loop order. Section[V]summarizes the obtained 
results. 



II. MODEL AND SCALING LAWS 

We consider a d-dimensional semi-infinite O(N) spin sys- 
tem which configuration is given by the iV-component clas- 
sical vector field s(r) satisfying the fixed-length constraint 
|s(r)| 2 = 1. The position vector r = (x, z) has a (d — 1)- 
dimensional component x parallel to the surface and a one- 
dimensional component z > which is perpendicular to the 
surface z = 0. It is convenient to introduce short notations 
for the volume integral over half space J v '■= J a dz J d l x 
and for the surface integral J s := J d d ~ 1 x. The large-scale 
behavior of the disordered spin system can be described by 
the effective Hamiltonian 



U [s] = Ho [s] + H suli [s] + H dis [s] 



(1) 



consisting of the sum of three terms which result from the 
semi-infinite bulk, surface and disorder in the bulk. The con- 
tribution from the the semi-infinite bulk can be expressed in 
the form of the well-known O(N) nonlinear sigma model: 



«o[s]= / 
Jv 



^(Vs(r)) 2 



h • s(r) 



(2) 



where h is the magnetic field in the bulk. The surface con- 
tribution to Hamiltonian can be written in its simplest form 
as IH: 



[s] 



hi • s(x), 



(3) 



where for simplicity we assume that the surface magnetic field 
hi has the same direction as the bulk field h. We consider a 
quite general type of bulk disorder such that its potential can 
be expanded in spin variables as follows 



EE 

/i=l ix-.-iu 



h^\(v) Sll (v)...s H Xv). (4) 
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The coefficients h) , (r) are assumed to be Gaussian ran- 

l \ . . . I f_ L \ > 

dom variables with zero mean and variances given by 



'•S.i.WCiM = ^Si^-Si^Sir - r'). (5) 

The first two coefficients have simple physical interpretation: 
h\ is a random magnetic field and h\j is a second-rank ran- 
dom anisotropy. The higher order coefficients h^) are higher 
order random anisotropics . As was shown in Ref. |2l] even 
if the system has only finite number of nonzero bare 
the RG transformations will generate an infinite set of higher- 
order anisotropies. However, the RG flow preserves the sym- 
metry with respect to rotation s — > — s. For instance, starting 
from the bare model with only a second-rank anisotropy only 
even-rank anisotropies will be generated by the RG flow. We 
will reserve the notation RA for the systems which possess 
this symmetry and the notation RF for the systems which do 
not. 

We employ the replica trick to average over disorder. Intro- 
ducing n replicas of the original system and averaging their 
joint partition function over disorder we obtain the replicated 
Hamiltonian 



i(Vs a (r)) 2 -h.s Q (r) 



1 U I 71 f 

2 T E ft(s«(r)-s 6 (r)) \ - £ / h x • s a (x), (6) 

a,b=l J a=l Js 



where we have defined the function IZ(z) = ^ 



The 



. _ v „, fJ „ 

properties of the original disordered system (flj can be ex- 
tracted in the limit n — > 0. According to the above definition 
of the RF and RA models, the function IZ(z) is arbitrary in 
the case of the RF model and even for the RA systems. 

Power counting shows that d\ c — 4 is the lower critical 
dimension of the model © B20I1 . Above the lower critical 
dimension the RF and RA systems undergo a paramagnetic- 
ferromagnetic transition. The scaling behavior at criticality 
is controlled by a zero temperature fixed point (FP) similar to 
the RFIM [45], reflecting the fact that disorder dominates over 
the thermal fluctuations. However, the temperature is danger- 
ously irrelevant. For instance, this results in violation of the 
usual hyperscaling relation and appearance of an additional 
universal exponent 9 that modifies the hyperscaling relation 
toJl: 



v(d-6) 



(7) 



where v and a are the correlation length and the specific heat 
exponents. One also expects a dramatic slowing down as the 
transition is approached with the characteristic relaxation time 
lnr ~ tT v , where t\ = \T — T c \/T c is the reduced temper- 
ature B46I1 . The magnetization in the bulk and on the surface 
vanish at the transition according to 



<r(ti)~tf, o-i(*i)~if, 



(8) 



where we have introduced the bulk and the surface magnetiza- 
tion exponents. At the critical point t\ = a small magnetic 



field in the bulk h induces the magnetization in the bulk and 
also on the surface according to 



a(h) ~ h 1 / 5 , cri(ft) ~h x ' Sl , 



(9) 



where we define the exponents i5 and Si . The surface magnetic 
field hi leads to the surface magnetization 



<7i(/ii) ~ h l 



l/<5i 



(10) 



Below the lower critical dimension d\ c a QLRO phase with 
zero magnetization can emerge. At criticality or in the QLRO 
phase, the correlation functions of the order parameter exhibit 
scaling behavior. Due to dangerous irrelevance of the tem- 
perature the connected and disconnected correlation functions 
scale with different exponents. We define the connected and 
disconnected correlation functions of the two local operators 
A and B as 



[A(r) ■ B(r')U := (A(v) ■ B{v')) ( A(v)) ■ ( B(r% 
[A(r) ■ B(r% is := (A(r)) ■ <B(r')> - (A(r)) ■ (B(r')). 

Here the angular brackets denote the thermal averaging and 
the bar stands for the disorder averaging. For instance, the 
connected and disconnected correlation functions of spins in 
the bulk scale independently as 



[s(r) • s(r')] con 
[s(r)- S (r')]dis 



1 



|r — r' \d-2+ri 1 

1 

| r _ r / 1 d— 4+7? ' 



(ID 

(12) 



Following the general scaling picture of the surface critical 
phenomena we introduce the surface exponents rj± and fj± 
which replace the bulk exponents 77 and fj in Eqs. ( fTTT ) and 
(fT2l) when one of the points r or r' belongs to the surface: 



[s(x, z) ■ s(x',0)] co 
[s(x,z) ■ s(x',0)]dii 



((x - x') 2 + z 2 ) (d - 2+v±)/[ 
1 

((x-x') 2 + ^ 2 ) ( ^ 4+,1)/2 ' 



:,(13) 



(14) 



We also define the surface exponents ryy and m which describe 
the connected and disconnected correlation function when the 
both points lie on the surface: 



[s(x) • s(x')] con 
[s(x) • s(x')] dis 



| x _ x /|d-2+r,n > 
1 

| x _ x ,|d-4+^ ■ 



(15) 



(16) 



Schwartz and Soffer ll47ll showed that the bulk exponents 
of the RF model obey the inequality 2rj > fj. The same ar- 
guments can be also applied to the surface correlation func- 
tions so that the surface exponents satisfy similar inequalities: 
277j_ > fj± and 2r]\\ > fj\\ . Note, that these inequality cannot 
be applied to the RA model where the coupling to disorder is 
bilinear. 
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III. FUNCTIONAL RENORMALIZATION GROUP 
A. Perturbation theory 

In the limit of low temperature and weak disorder the con- 
figuration of the system is fluctuating around the completely 
ordered state in which all replicas of all spins align along the 
same direction which is parallel to h and h x . It is convenient 
to split the order parameter s a = (a a , 7r a ) into the (N — 1)- 
component vector n a which is perpendicular to this direction 
and the component a a = -Jl — tt 2 parallel to this direction. 
Then the effective action of the system can be written as 

1 a , fa, ' V7T a ) 2 



1 n 



fa 

2(1 - ttD 



ha a 



7T5 + cr a (Th) . (17) 



In general one has to add to the action ( fTvT > the terms like 
5 d (0) J v ln(l — tt 2 ) generated by the Jacobian of the trans- 
formation from s a to 7r a . However, in what follows we 
will use the dimensional regularization scheme ll48ll in which 
8 d (0) = so that we ignore these terms in action (TTTb from 
the beginning. 

Let us denote averaging with the action (FT/l by double an- 
gular brackets and introduce the following correlation func- 
tions 



where the connected correlator corresponds to a single replica 
and the disconnected one to two different replicas a ^ b. To 
compute the correlation functions at the surface like [71- (r) ■ 
7r(x')] con or [7r(x) • 7r(x')]di s one has to replace G^ 2,0 * 1 by 
G* 1 ' 1 ) and G^°' 2 \ respectively. 

Expanding the effective action ( fT7b in small tt we will treat 
the quadratic part as a free action and the rest infinite series 
as interaction vertices (see Appendix [A}. Then the correla- 
tion functions ( fT~8T > can be expressed in terms of Feynman 
diagrams which give the low temperature and small disor- 
der expansion. In practical calculations it is convenient to 
perform the Fourier transform with respect to x: 7r(q, z) 
f d d - 1 xTv(x,z)e- lcl - x and define / := J cF^VO") 
The quadratic terms give the free propagator 



\d-l 



1 



e ~q\z-z'\ + 



h _ 



q(z+z') 



hi 



(23) 



where we have introduced the short notation q := (q 2 + h) 1 ^ 2 . 
The free propagator ((23) satisfies the boundary conditions 



[d z -/ii]G(°>(x,z,x',z') 



0. 



z=0 



(24) 



The free surface corresponds to the limit hi — > in which 
Eq. d23l becomes the Neumann propagator consisting of the 
bulk part and the image part. In what follows we will use 
the Neumann propagator as the bare one and treat the terms 
proportional to hi as soft insertions $44 



where L points r = (ri, are off surface and K points 

x = (xi, ...,Xk) are siting on the surface. In Eq. (fT8T > we 
have used a short notation a = (cei, «l) where each a u 
stands for the component number i v and the replica number 
a„. The similar holds for /3, The correlation functions (fT~8t can 
be computed using the following generating functional [49] 

JF[J,Ji] =ln J 2? 7re - sw+ ^ J ( r ) 7r ( r )+/s J i( x ) 7r ( x ), (19) 

where we assume that the source J(r) vanishes at the surface. 
Differentiating with respect to the sources we obtain 

«x>=n^n^| ,<*> 

v ~ 1 t 1 - 1 J=Ji=0 

where for the sake of brevity we have suppressed all tensorial 
indices. Using correlation functions (fT~8t one can compute 
the connected and disconnected functions defined in Eqs. ( fTTI l 
and ( fT2l . However, since we are interested only in the scaling 
behavior it is more convenient to consider the similar corre- 
lation functions not for s but for tt fields. For example, the 
correlation functions at two off surface points read 



N-l 



[7r(r).7r(r')] c 



lim VGg$ a (r,r'), (21) 



[7r(r) • 7r(r')] dis = lim 



+0 ^— ' 

i=l 

JV— 1 



^ 

1 



^(2,0) / /s 



(22) 



B. FRG equations and critical exponents 

The correlation functions ([T8T l calculated perturbatively 
in small disorder and temperature suffer of the UV diver- 
gences. To avoid mixture with IR singularities in the O(N)- 
noninvariant correlation functions it is convenient to keep 
h ^ 0. The UV divergences can be converted into poles in 
e = d — 4 using dimensional regularization. To renormalize 
the theory one has to absorb these poles into finite number of 
Z-factors. However, all the Taylor coefficients of the dis- 
order correlator TZ(<f>) turn out to be relevant operators so that 
one has to introduce renormalization of the whole function. To 
simplify calculation of the disorder renormalization one can 
use the background field method rf25tl . Using the Legendre 
transform of the generating functional ([19} from the sources 
J to the background fields LT one derives the effective action 
r [LT] which is the generating functional of the one-particle ir- 
reducible vertices. The two-replica part of the effective action 
gives the renormalization of the disorder. Since the scaling 
behavior is controlled by a zero temperature FP we will dis- 
regard all terms involving more than two replicas which are 
suppressed in the limit T — >• 0. The renormalization of the dis- 
order simplifies by changing variables: lZ(<fi) — R(z) where 
z = coscf), for instance, TZ'(1) = —R"(0). In terms of the 
variable <fi, the function R{<j>) becomes periodic with the pe- 
riod 2ir in the RF case and with the period tt in the RA case. 
The relation between the renormalized and the bare correla- 
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tion functions reads 

G^ (r; T, h, h lt R, /i) = Z^ L+K ^ 2 Z; K/2 

where circles denote the bare quantities and \i is an arbitrary 
momentum scale. UV divergences are absorbed into Z-factors 
according to 

* = Zy 2 7r, tt\ s = {Z^Zx) 1 / 2 ir\ s , (26) 

h = n 2 Z T Z- 1/2 h, h 1 =nZ T (Z ir Z 1 )- 1 ' i hi, (27) 
T = /i 2 - d Z T T, R = n 4 - d K7 1 Z R [R], (28) 



where {2ir) d K d = 2n d/2 /T(d/2) is the surface area of a d- 
dimensional unit sphere and T(x) is the Euler gamma func- 
tion. Zr[R] in Eq. d28b is a functional acting on the renormal- 
ized disorder correlator R(<j>) which has the following loop 
expansion: 



Z R [R] = R + S W (R,R) + S (2) (R 1 R,R) 



(29) 



where S^(R,R) is bilinear in R and proportional to 1/e, 
while 5^ [R, R, R) is cubic in R and contains terms of or- 
der 1/e and 1/e 2 . According to Eq. d26l i the surface field tt\ s 
renormalizes differently from the field tv in the bulk. The new 
factor Z\ serves to cancel the additional UV divergences in 
Feynman diagrams arising from the image part of the Neu- 
mann propagator Gq (z,z') for z' — > 0. The renormal- 
ized theory is not unique and depends on the scale /i. Using 
this fact we will derive the functional renormalization group 
(FRG) equation. 

We now consider how the scaling behavior can be extracted 
from the renormalized theory. Using independence of the bare 
theory on the momentum scale fi one can derive the flow equa- 
tions for the renormalized correlation functions differentiating 
the both sides of Eq. d25l l with respect to fj, at fixed bare quan- 
tities. One finds that the renormalized correlation functions 
satisfy the following FRG equation 



M 3 M + (d - 2 - Cr)Td T - C h hd h - Ch.hxdh, + 



q{l,k) = Q) (30) 



where the integral in the last line is taken over a period, i.e., 
(0, 7r) for RA and (0, 2tt) for RF models and we have intro- 
duced the scaling functions: 

C, = fj,d M ]nZi\ , (i = T,%,l), (31) 

C h = 2 + Ct - CrA (32) 

On =1 + Ct-(Cx + Ci)/2, (33) 

P[R] = -Hd li R(<l>)\ Q . (34) 

Here the zero indicates that the derivatives are taken at fixed 
bare quantities. Flow equations similar to Eq. (l30t hold also 
for the correlation functions in which some or all the fields 
7r Q (r) are replaced by a a (r) and for other observables, e.g., 
the correlation length and the magnetization [48]. 



The long-distance physics can be obtained from the solu- 
tion of the FRG equation fl30l ) in the limit of /j, — > 0. The 
renormalized disorder correlator and the temperature flow ac- 
cording to 



-fidfj, In T = 2-d 



(35) 
(36) 



The scaling behavior is controlled by a zero temperature FP 
/3[R*} = with R* of order e and T* = 0. Indeed, according 
to Eq. d36l >. the temperature is irrelevant, i.e. it flows to in 
the limit /i — > for d > 2 and for sufficiently small £x = 
O(R). Although one expects that (t is small in the vicinity of 
the FR one has to take caution whether the zero temperature 
FP survives in three dimensions where (r ~ e is negative 
l26Tl . The stability of the FP can be checked by computing 
the eigenvalues of the disorder flow equation (|35l l linearized 
about the FP solution: R{<j>) = R*(<jj) + J2i k%(4>)- Since 
one expects that for d > 4 (e > 0) the FP R* (tp) describes the 
paramagnetic-ferromagnetic transition it has to be unstable in 
a single direction ^i(<j>) with eigenvalue Ai > 0: /3[R* + 
ti^i] = Aiti^i + 0(t 2 ). In vicinity of the zero temperature 
FP which controls the paramagnetic-ferromagnetic transition, 
the FRG equation for the correlation length £ can be written 
as 



d_ 



£(/*,ti) = 0. 



(37) 



Dimensional analysis implies that t\) = [i 1 <^(ti). This 
reduces Eq. d37| i to an ordinary differential equation (ODE) 

which solution is £ ~ [i~ 1 t l 1 ^ Al . The latter describes di- 
vergence of the correlation length on the critical line at zero 
temperature when the strength of disorder approaches the crit- 
ical value B45I1 . Assuming that along the transition line at fi- 
nite temperature t\ ~ T — T c we find that the positive eigen- 
value Ai gives the critical exponent of the correlation length 
v = l/Ai. For d < 4 (e < 0) the FP becomes stable and 
describes a QLRO phase. The fluctuations exhibit power-law 
correlations in the whole QLRO phase so that the correlation 
length £ is always infinite down to the lower critical dimension 
of the QLRO - paramagnetic transition. 

Let us consider the solution of Eq. (l30T > for the connected 
two-point correlation functions. The dangerous irrelevance 
of the temperature manifests itself in the fact that the con- 
nected (bulk or surface) two point functions are proportional 
to T in the low temperature limit. This is explicitly shown 
in Appendix lAl for the connected correlation function G^ 1 ' 1 ). 
Hence, setting h — hi = and R = R* we can rewrite 
Eq. (O as 



^ + i(L+if)c; + |cr+< 



G{ 



L.K) 



0, (38) 



where the star denotes that the function is computed at the FP. 
In Eq. ( 1381 ) we have defined the exponent 



d-2-C 



T' 



(39) 



which describes the flow of the temperature (136) in the vicin- 
ity of the FP and which has been introduced ad hoc in the 



6 



modified hyperscaling relation Using the method of char- 
acteristics and dimensional analysis one can write the solution 
of Eq. ( T38t in the form 

Gi L ok K \rb- R*) = b-^ L+K ^ +K ^' 2 +^f c {r-R*). (40) 

Considering the connected two point functions d40l i with (L = 

2, K = 0), (L = 1, K = 1), and (L = 0,K = 2) we derive 
the critical exponents: 



v\\ = £ + 



(41) 
(42) 
(43) 



We next turn to the disconnected two-point correlation 
functions. At variance with the connected correlation func- 
tions they are not proportional to the temperature. Thus, at 
h = hi = T = they satisfy the same Eq. ( f38l > but with- 
out the term 8 in the brackets. The solution of the latter FRG 
equation is given by 

G ( jQ K \rb;R*) = b-^ L+K ^ +K ^^ f d (r; R*). (44) 

Repeating analysis we did for the connected functions we ar- 
rive at 

f} = 4-d + £ = 2 + V -6, (45) 
fjx=4-d + C + &72 = 2 + r a -0 1 (46) 
rj||=4-(f+C + Ci =2 + 1711-0. (47) 

Note that the exponents d4Tb-(f43l> and (|45l>-(f47b are related by 



2t7_l = V + i7|| , 2rj± = 7) + r]\\ 



(48) 



Finally we study the profile of the spontaneous magnetiza- 
tion below and at the paramagnetic-ferromagnetic transition 
for d > d\ c . The magnetization as a function of the distance 
to the surface z, the reduced temperature t\, and the bulk and 
surface magnetic fields h and hi satisfies the following flow 
equation 



d 

oti 



a(z,tx,h,hi) = 0. (49) 



Here j = and z > corresponds to the bulk magnetiza- 
tion a while j = 1 and z = gives the surface magnetiza- 
tion <j\ . The solution of Eq. d49l ) can be written as 

£r(z,ti,Mi) = &-(4«+* c ?> 

xcr(z^ 1 , tifc Al , /i& c '* , /ii6 c; i ). (50) 

We first consider the profile for h = hi = 0. The solution 
( |50l interpolates between the surface magnetization <Ji(ti) ~ 
t (C+CD/(2Ai) atz ~ o and the bulk magnetization cr (ii,z) ~ 
for z 3> £■ Reexpressing the latter in terms of ^,77, 



C/(2Ai) 



f 

and 7711 we obtain that the bulk and the surface magnetization 
exponents defined in Eq. ((H) are given by 

= ±u(d-A + fj), pi = ^v(d-4 + fj l{ ). (51) 




<8 JtWWWwA • 

(b) 



(c) 



FIG. 1 . One-loop diagrams contributing to the connected two-point 
function G^'^-l a (z, p). The solid lines stand for the Neumann prop- 
agator ([23}. The wavy and dashed lines are vertices defined in 
Eqs. JAH >- fA3l >. The crossed circles denote the points on the surface. 



At the critical point ti = and finite external fields we find 
that o-(/i) - ft,C/( 2 a) in the bulk and ai (h) ~ fo(C+CD/(2«) 

or <Ti(hi) ~ hi' + ^^ at the surface. Thus, the expo- 
nents 6, Si, and 6u defined in Eqs. (0 and ( fTOb satisfy the 
following scaling relations: 

6-1 v 81- P/fii v 811 - 1 



2-17 /3' 



2-17 



Pi 



1 



(52) 



IV. THE SURFACE EXPONENTS TO ONE-LOOP ORDER 

We now renormalize the both semi-infinite RE and RA 
models to one-loop order and explicitly calculate the surface 
critical exponents to first order in e = d — 4. The factors 
Zn, Zt and Zr[R] defined in Eqs. (|26T>-(|29T> are the same that 
appear in the case of the infinite systems. They have been cal- 
culated in several works up to two-loop order 1I2TI I26I - I29I . To 
one-loop order they read 

Z« = 1-(N-1^ 



Z T = 1 - (N - 2) 



£ 

R"{0) 



1 



-0(R 2 ), (53) 
hO(.R 2 ), (54) 
R"(0)R"(<t>) 
(N - 2)iR"(0)[2R(<p) + R'{(j)) cot cj>] 

1 [R'(cf>)} 2 \. (55) 



sS^(R ) R) = -R"(<f>) 



2 sin 2 , 

The new factor Z\ which eliminates the poles resulting from 
the presence of the surface can be determined from the renor- 
malization of the two point function (p, z; h, T, R). The 
one-loop diagrams contributing to this function are shown 
in Fig. Q] The corresponding integrals are computed in Ap- 
pendix [A] and give 

GW(p,z;h,f,R) = f^- \l-^R"(0) 



P 



4e 



(iV-3) + — I + 2(tf + l) 
\P P 



0(R 2 ) } , (56) 



where p = (p 2 + h 2 ) 1 / 2 . The factor Zi can be found from the 
renormalization condition 

Z^Z^^G^Hp.z-Xf ,R) = finite for e 0,(57) 



7 



where the bare h, f, R are replaced by the renormalized h, T 
and R according to Eqs. (|26]l-(|28T>. We obtain 

Zi = l-(N- + 0(R 2 ). (58) 

e 

Thus, to one loop order we have Z\ = + 0(R 2 ). Using 
Eqs. (|3~TT > and (l34l i we calculate the scaling functions 

(t = -{N-2)R"(0) + O(R 2 ), (59) 
C w = Ci = -(iV-l)ii"(0) + O(i? 2 ), (60) 



and the beta function 



1 



P[R] - -eR(cj>) + ^R"(4>f - i?"(0)i?» 



-(N - 2)|i?"(O)[2i?(0) + R'{4>) cot . 
' [i?'(0)] 2 )+O(i? 2 ) 



2siri 



(61) 



to one-loop order. Solution of the FP equation f3[R*} = 
with the beta function (|6Tb has been analyzed for different 
values of N and different sign of e in Refs. I26l - l29ll . We first 
assume for granted that the flow has a FP R* ((f)) which is a ir- 
periodic function for the RA model and a 2ir -periodic function 
for the RF model. Then, the surface critical exponents can be 
computed to one loop using Eqs. d4TT)-(l43b and d45T>-(l47T> that 
yields 



n 



-R*"(Q), 
N + l 



V 



-s-(N-l)R*"(0), (62) 



V± = —R*"(0), fj ± = -e- -(N - 1)R*"(0), (63) 



m = -NR*"(0), 



fju =-e- 2(N - 1)R*"(0). (64) 



by the 



The other surface exponents are related to 
scaling relations dBTl and (l52l . 

Before we explicitly calculate the surface exponents for the 
semi-infinite RF and RA models let us remind how the FRG 
allows one to overcome the DR problem. The incorrect DR 
prediction results from the assumption that the flow equation 
(l35T l with the beta function doTt has a FP which is an analytic 
function. Indeed, in this case one can obtain a closed flow 
equation for the R"(0): 



^R"{0) = -sR"(0) — (N — 2)R"(0) 2 



(65) 



Equation d65l > has a nontrivial FP solution i?*"(0) = 
—e/(N — 2) with the eigenvalue Ai = e. This FP is un- 
stable for e > as one expects for a FP corresponding to the 
transition and gives the DR exponents: j/ DR ) = 1/e and 



(DR.) 



(DR) 



7Y-2' 
N+l 



2(N 
N 



(DR) _ (DR) _ . 

^11 ~N-2- 



(66) 
(67) 

(68) 
(69) 



The one-loop DR exponents for the magnetization read 







(DR) 



N- 1 



^(DR) = A 



1 



(70) 



2(iV-2)' " 1 TV -2' 

For e < the FP is stable but the rj critical exponents become 
negative, and hence, unphysical. 

More accurate analysis of the RG flow shows that R"'(0) 
diverges at a finite scale /i. Thus, no analytic FP can exist 
and one has to look for a non-analytic FP with R*"'(0 + ) ^ 
which would violate the DR predictions. This requires so- 
lution of the boundary-value problem for the nonlinear ODE 
with periodic boundary conditions, which depend on the uni- 
versality class. We assume that the small <f> expansion of the 
FP solution R*(4>) has the following form 



R*(4>) = a + a 2 <t> 2 + a 3 |0| 3 + a 4 </> 4 + a 5 |0| 5 + 



(71) 



meaning that R*"(<fi) has a cusp at the origin with R*"'(Q + ) ^ 
0. Substituting ansatz ( TTTb into the FP equation we find that 
the first coefficients are given by 

2a 2 (N-l) _ R*"(0) 



a 



«3 



4(N -2)a 2 + e' 



« 2 



-sign(e) 1 



' 2ea 2 + Aa 2 (N - 2) 



9(^ + 2) 



(72) 



(73) 



The value of R*"(0) as well as the sign of 03 are constrained 
by the boundary conditions. R*"(0) can be determined us- 
ing the shooting method to fulfill the appropriate periodicity 
requirement. 



A. Random field O(N) model 

1. Paramagnetic-ferromagnetic transition for d > 4 (e > 0) 

The RF model is described by R((f>) which is a 27r-periodic 
function. Numerical solution of the FP equation shows that 
for d > 4 a 2ir -periodic solution of the form (l7ll-(l73l exists 
only for N > N c — 2.834 74. It has R*"(0) < and it dis- 
appears when — > 7Y+. This cuspy FP is once unstable with 
the positive eigenvalue Ai = e. Thus, the correlation length 
exponent v = 1/e + (e°) coincides with the DR prediction to 
one-loop order. Remarkably, the non-zero R*"'(0 + ) vanishes 
for N > N* = 18 + 0(e). The non-analyticity becomes 
weaker as N increases and starts with i?*( 2 P( A ')+ 1 )(0 + ) 
where p ~ N |28, 29, 52]. Weaker non-analyticity results in 
restoring the DR critical exponents for N > N* . The crit- 
ical exponents rji and fji computed using Eqs. d62l-(l6*4l) as 
functions of N are shown in Fig. [2] With increasing N they 
monotonically decay approaching the DR values at N = N* 
and satisfying the inequalities: rj < fj < rj± < fj± < m < fjn. 
The bulk and surface magnetization exponent j3 and /3\ calcu- 
lated for different N are shown in inset of Fig. [2] To one-loop 
order they obey relation j3i — 2(3. Up to now the both mag- 
netization exponents have been studied only for the 3D RFIM 
where numerical simulations give /3 = 0.0017 ± 0.005 fBUl 
and fa = 0.23 ± 0.03 JH. Thus, the ratio fo/p for the 
RF O(N) systems in d > 4 is much smaller than for the 3D 
RFIM. 
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FIG. 2. (Color online) The critical exponents rji and fji (divided by 
e), which describe the paramagnetic-ferromagnetic transition of the 
RF model above the lower critical dimension, as functions of TV for 
TV > N c . Inset: The corresponding bulk magnetization exponent ft 
and the surface magnetization exponent /3i as functions of TV. 




N 



FIG. 3. (Color online) The critical exponents rji and fji (divided by 
je|), which describe the power-law decay of correlations in the QLRO 
phase of the RF model below the lower critical dimension, as func- 
tions of TV for TV < N c . 



2. Quasi-long-range order for d < 4 (e < 0) 

Below the lower critical dimension the flow equation for the 
disorder correlator has an attractive 2vT-periodic FP solution of 
the form (1711 -(1731. This cuspy FP appears only for 2 < N < 
N c where it controls the scaling behavior of spin fluctuations 
in the QLRO phase. The corresponding exponents rji and fji 
as functions of N are shown in Fig. [3] In the case N = 2 the 
FP equation admits for an explicit non-analytic 0o -periodic 
solution given by 



iT(0) 



72 



1 

36 



Using Eqs. (f62t-(l64b one obtains 



V = 

v\\ = 



36 1 
24 1 
18 1 



V 



2 / i \ 2 

00 



3G 



f )\\ = 



1 



0| 

24 

00 

18 



(74) 

(75) 
(76) 
(77) 



with O = 27r for the RF system. 

The semi-infinite RF 0(2) model can be mapped onto a 
semi-infinite periodic disordered elastic system with a free 
surface. There is one to one correspondence between the 
Bragg glass phase of the elastic system and the QLRO phase 
of the studied spin model. The power-law decay of the spin 
correlations in the QLRO phase corresponds to the logarith- 
mic growth of the displacements in the disordered elastic sys- 
tem. Moreover, the exponents r\, rj± and m provide the uni- 
versal amplitudes of the logarithmic growth of the displace- 
ments in the bulk, at the surface and along the surface, respec- 
tively. For a 0o -periodic elastic system with a free surface 
these amplitudes are given by Eqs. (|75ll-(|75ll. In particular, we 
find that the logarithmic growth of the displacements along the 



surface is twice larger than the logarithmic growth in the bulk. 
In the case when only one point is on the surface the growth is 
enhanced by 50%. The presence of a free surface can be con- 
sidered as an extended defect of a special kind. The influence 
of potential-like extended defects on the Bragg-glass has been 
recently studied in Refs. HIM]. 



B. Random anisotropy O(N) model 

1. Paramagnetic-ferromagnetic transition for d > 4 (e > 0) 

The FP equation has a cuspy 7r-periodic solution of the 
form (f7lTi-(l73l which is singly unstable giving the correla- 
tion length exponent v = 1/e + (e°). It exists for any 
N > N c = 9.4412 with a non-zero R*'"(0+). Therefore, 
at variance with the RF case in the RA model the DR breaks 
down for all values N > N c , i.e., N* = oo HI]. The N- 
dependence of the critical exponents r\; L and fji is shown in 
Fig. |4] For large N one can find the asymptotic behavior of 
the FP solution Q Si HI El- Following Ref. JH we look 
for the 7r-periodic solution of the FP equation @\R\ =0 with 
the beta function doTT l of the form 



R*'{4>) = --Seam 



IT -2(f) 



(2x(0) - l)G(x). (78) 



Here we have introduced a small parameter 5 = 1/(N— 2) and 
defined variable x(4>) = cos(^p^). Substituting ansatz ( F78b 
into the FP equation and expanding the function G{x) in small 
5 one finds that the coefficients are polynomials in x: 

G{x) = 1 + ^(95 - AAx - l%x 2 )5 - —(11737 - 5040a; 
9 81 



-3624a; 2 - 3104x 3 - 768x 4 )c5 2 



10935 



(103378933 



-45854072a; - 23128624a;^ - 16172328a;' :i 

-9791216a; 4 - 4642048a; 5 - 901120x 6 )^ 3 + 0{8 A ). (79) 
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FIG. 4. (Color online) The critical exponents r\i and fji (divided by 
e), which describe the paramagnetic-ferromagnetic transition of the 
RA model above the lower critical dimension, as functions of N for 
N > N c . Inset: The corresponding bulk magnetization exponent /3 
and the surface magnetization exponent /3i as functions of N. 
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FIG. 5. (Color online) The critical exponents r\i and fji (divided by 
|e|), which describe the power-law decay of correlations in the QLRO 
phase of the RA model below the lower critical dimension, as func- 
tions of N for N < N c . 



This implies that ||53| 

Substituting the solution (l80l into Eqs. do*2l-(l64l we find the 
correlation function exponents to leading order in 1 /N as: 



3e / 52 \ _ e / 49 



(81) 



3e / 55 \ 5e / 147 

^ = T 1 + 5aF + - »^ = T 1+ saF + - - (82) 



r )\\ 



4 
3c 



3N 
52 



, 1 + — + •■■ , m = 2e 1 , 



oN 
49 



H 4 V 3JV 



H 2 I 37V 



(83) 
(84) 



where in the last line are the bulk and the surface magnetiza- 
tion exponents. 



models have the lower critical dimension d\ c = 4. Above 
die they undergo a paramagnetic-ferromagnetic transition for 
N > N c , while below d\ c and for N < N c they exhibit a 
QLRO phase with zero magnetization and power-law correla- 
tion of spins. Here the critical value N c = 2.835 for the RF 
models and N c — 9.441 for the RA. Using FRG we studied 
the surface scaling behavior of these models at criticality as 
well as in the QLRO phase, and calculate the corresponding 
surface exponents to lowest order in e = d—4. We have found 
that the DR prediction for the surface scaling is broken similar 
to that happens in the bulk. We have shown that the connected 
and disconnected correlation functions scale differently also 
at the surface and derived the scaling relations between dif- 
ferent surface exponents. The surface exponents obtained for 
the 3D RF 0(2) can be used to describe the growth of dis- 
placements near a free surface in semi-infinite periodic elastic 
systems in disordered media. The surface scaling we obtained 
for the Heisenberg {N = 3) RA model can be relevant for the 
behavior of amorphous magnets iTToL l55ll . 



2. Quasi-long-range order for d<4(e<0) 



For 2 < N < N c the flow equation has a stable 7r-periodic 
FP solution of the form d7ll-(f73l which controls the scaling 
behavior of spin fluctuations in the QLRO phase of the RA 
model for d < 4. The correlation function exponents rji and fji 
computed for different N are shown in Fig. [5] For N = 2 the 
FP equation has an explicit non-analytic 7r-periodic solution 
given by Eq. (l74l i with <po = it. The critical exponents of the 
RA 0(2) model are given by Eqs. (|75]l-(f77]i with </> = n. 
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V. SUMMARY 



Appendix A: One-loop diagrams contributing to 



In the present work, we have investigated the RF and RA 
semi-infinite O(N) models with a free surface. The both 



In this appendix we calculate the correlation function 
G^ia-i a t0 one-loop order. Expanding action (fTTt in small 
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7r a we find that the only vertices we need are 





ia, 




ia 


ia 




Ja 

z 


-> 

z 


= -4-^"(o), 

2T 2 V " 


J a 

> 


jb 


= — i-£"(o). 

8T 2 



g 2 + a zi a z 



,(A1) 
(A2) 
(A3) 



The one loop diagrams contributing to the correlation function 

a ^ shown in Fig. Q] The solid line corresponds to 
the Neumann propagator (|23l with hi = and the wavy and 
dashed lines to vertices (IAlb -( |A3l . The first diagram gives 



loop integral reads 
1 



q{p + q) 2 



q d ~ 2 dq 



(q 2 + h 2 ) 3 / 2 



O(e 



IK, 



- + 0(e°). 



(A5) 



The second and third diagrams yield 
1 



oc oo oo 



(b) = j^R"{0) J J dz l J dz 2 J dz 3 S(z 2 - zi) 
q 

(p + q) 2 + d Zl d Z2 +h] GW(0, Zi )GW(z 2 ,z) 

x G^(z 1 ,z 3 )G^(z 3 ,z 2 ) = -l-ffr'We-** 



zh h 

— + - 



2pz-l\l 2 - 2I 3 



finite, (A6) 



oo oo oo 



(a) = —^1r"(0) J J d Zl J dz 2 J dz 3 S(z 2 - z x ) 
9 

d Zl d Z2 +h] G^{Q, Zl )G^{ Zl ,z) \gW(z 2 ,z 3 ) 



^—^fR"(0)e-f z [ — + % + 2 ) I 2 + finite,(A4) 
lop \ P P ] 



where we have used p — (p 2 + ft, 2 ) 1 / 2 and omitted the terms 
finite in the limit e — > 0. The logarithmically divergent one- 



(c) = -±R"(0) J j dz 1 G^(0,z 1 )G^(z 1 ,z 1 ) 
9 o 

xG^( Zl ,z) = -^TR"(0)e- pz 
x [(2pz + 5)I 2 + 2I 3 ] + finite, (A7) 
where we have defined the algebraically divergent integral 

3p + q + (2p + q)pz 



h(p,z) 



p 2 (p + q) 2 



(A8) 



Summing up the three diagrams we find that the algebraically 
divergent integral (IA8b cancels and we obtain Eq. 
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